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Inspirals of neutron star-neutron star binaries are a promising source of gravitational waves for 
gravitational wave detectors like LIGO. During the inspiral, the tidal gravitational field of one of 
the stars can resonantly excite internal modes of the other star, resulting in a phase shift in the 
gravitational wave signal. We compute using a Fisher-matrix analysis how large the phase shift 
must be in order to be detectable. For a 1AMq,1AMq binary the result is ~ 8.1,2.9 and 1.8 
radians, for resonant frequencies of 16, 32 and 64 Hz. The measurement accuracies of the other 
binary parameters are degraded by inclusion of the mode resonance effect. 
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I. INTRODUCTION AND SUMMARY 

In neutron star - neutron star (NS-NS) binaries, each 
star exerts a tidal gravitational force on its compan- 
ion as the stars inspiral. This force can resonantly 
drive internal modes of oscillation of the companion 
0, H 0, 1, H, 0, B !, El, which in turn can alter the 
emitted gravitational wave signal. The mode driving can 
be primarily dissipative, or adiabatic, or resonant; see 
for example the discussion in Ref. [9(. The modifica- 
tion to the emitted gravitational wave signal can poten- 
tially yield information about the internal structure of 
the stars, or impede the matched-filtering based detec- 
tion of the events. 

In this paper we will focus attention on resonant mode 
driving. To a good approximation, the effect of the mode 
resonance on the gravitational wave signal can be de- 
scribed by two parameters, the gravitational-wave fre- 
quency fo (twice the orbital frequency) at which reso- 
nance occurs, and a constant phase shift parameter A$. 
If we denote by $o(/) the phase of the Fourier transform 
of the gravitational waveform neglecting mode driving, 
then the phase $(/) including mode driving can be writ- 
ten as H 



*(/) 



*o(/) + (1 - ///o)A$ f - / » A/ r , 
<f>o(/) / - Jo » A/ r , 



(1.1) 

Here we have chosen the parameters of the unperturbed 
and perturbed waveforms so that the waveforms coin- 
cide after the resonance. Also A/ rcs is the bandwidth of 
the resonance, which is sufficiently small 1 that we can 
neglect it for data-analysis purposes. Thus, the phase 
perturbation due to the resonance is a linear function 
of frequency which vanishes after the resonance, whose 
maximum value is ~ A$. 

There are examples known of situations for which the 



resonant phase shift A<f> is large compared to unity, which 
implies the effect should be easily detectable in the grav- 
itational wave signal. However all such cases requires the 
spin frequency of the star to be larger than is thought 
to be likely for most NS-NS inspirals. Resonant excita- 
tion requires that the mode frequency be small compared 
with the natural frequency oj ~ y/MjB? of the star, 
where M is the stellar mass and R the stellar radius. 
One class of modes with suitably small frequencies are 
g-modes; however the overlap integrals for these modes 
are so small that A$ is small compared to unity @, [H, 0] • 
Another class are the / and p-modes of rapidly rotating 
stars, where the inertial-frame frequency ui m can b e much 
smaller than the corotating-frame frequency ~ -J AI/R 3 . 
Ho and Lai [5[ showed that the phase shifts due to these 
modes could be large compared to unity. However, the re- 
quired NS spin frequencies are several hundred Hz, which 
is thought to be unlikely in inspiralling NS-NS bina- 
ries. A third class of modes are Rossby modes (r-modes) 
[TTL [l2j , for which the restoring force is dominated by the 
Coriolis force. For these modes the mode frequency u m 
is of order the spin frequency of the star, and thus can be 
suitably small [10 Hz < u m /{2n) < 100 Hz]. Ho and Lai 
|5| computed the Newtonian driving of these modes, and 
showed that the phase shift is small compared to unity. 
Ref. @ showed that a larger phase shift is produced by 
post-Newtonian, gravitomagnetic driving, but that once 
again A<i> > 1 is only possible for spin frequencies of 
order several hundred Hz. 

In this paper we compute how large the phase shift 
A$ needs to be in order to be detectable. We confirm 
the prevailing expectation that the detectability criterion 
is A$ > 1. More specifically, we consider a simplified 
model of the gravitational wave signal [l[ which depends 
only on 5 parameters: the masses Mi and Mi of the 
two stars, the time t c of coalescence, the orbital phase 
4> c at coalescence, and an overall amplitude parameter A 
2 . Table HI from Ref. [l| , shows the accuracy with which 
these parameters could be measured in the absence of 



1 From Eqs. (1.5) and (3.2) of Ref. @|, the bandwidth is 

M \ 5/(i ( fo ^ 11/6 



A/ r< 



■ 0.1 Hz 



,1.2 M & J VlOHz 
where M is the chirp mass of the binary. 



2 For simplicity we neglect the influence of the spins of the NSs on 
the parameter extraction. Including spin will degrade somewhat 
the measurement accuracies. 
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TABLE I: The rms errors for signal parameters, neglecting 
spin effects and in the absence of resonances, assuming the 
advanced LIGO noise spectrum [lj], and for a signal-to-noise 
ratio of 10. The units of the stellar masses Mi and M2 are 
solar masses, while At c is in milliseconds. 



Mi 


M 2 


At c 


A0 C AM/M Afi/n 


2.0 


1.0 


0.72 


1.31 


0.004% 0.39% 


1.4 


1.4 


0.71 


1.28 


0.004% 0.41% 


10 


1.4 


1.01 


1.63 


0.020% 0.54% 


15 


5.0 


1.44 


2.02 


0.113% 1.50% 


10 


10 


1.43 


1.98 


0.160% 1.90% 



TABLE II: The rms errors for signal parameters, as in Ta- 
ble HI but including the effects of resonances. The resonant 
frequency fo is in Hz. The quantity A(A$) is the rms mea- 
surement error in the phase shift parameter A"l>. 



Mi 


M 2 


fo 


At c 


A(j} c 


AM/M A/z//i 


A/o 


A(A$) 


2.0 


1.0 


16 


0.81 


1.68 


0.007% 


0.56% 


24 


8.1 


1.4 


1.4 


16 


0.80 


1.65 


0.008% 


0.58% 


24 


8.1 


10 


1.4 


16 


1.20 


2.22 


0.040% 


0.82% 


24 


8.6 


15 


5.0 


16 


1.82 


2.89 


0.232% 


2.44% 


24 


9.2 


10 


10 


16 


1.81 


2.85 


0.318% 


3.14% 


24 


9.2 


2.0 


1.0 


32 


0.76 


1.46 


0.007% 


0.47% 


17 


2.9 


1.4 


1.4 


32 


0.75 


1.43 


0.007% 


0.48% 


17 


2.9 


10 


1.4 


32 


1.10 


1.86 


0.033% 


0.64% 


18 


2.9 


15 


5.0 


32 


1.64 


2.36 


0.165% 


1.82% 


18 


2.9 


10 


10 


32 


1.61 


2.32 


0.224% 


2.33% 


18 


2.9 


2.0 


1.0 


64 


0.97 


2.44 


0.005% 


0.85% 


23 


1.8 


1.4 


1.4 


64 


0.95 


2.39 


0.005% 


0.88% 


23 


1.8 


10 


1.4 


64 


1.59 


3.47 


0.030% 


1.31% 


23 


2.1 


15 


5.0 


64 


2.48 


4.55 


0.202% 


3.87% 


23 


2.5 


10 


10 


64 


2.45 


4.45 


0.280% 


4.93% 


23 
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mode excitation, for an event with signal-to-noise ratio 
of 10, and assuming the advanced LIGO noise spectrum. 

Next, we enlarge the signal parameter space to in- 
clude the two parameters fo, the gravitational wave fre- 
quency at resonance, and A<1>, the phase shift parameter. 
We compute measurement accuracies for binaries with 
masses (2M ,1M ), (1.4M Q , 1.4M Q ), (1OM , 1.4M©), 
(15M©, 5M©), (10Mq,10M q ) and resonant frequencies 
fo = 16 Hz, 32 Hz and 64 Hz with A$=l. The results 
are given in Table [Til The entries in the last column of 
this table are the minimum values of A$ necessary for 
detectability of the resonance effect (see Sec.[TlJ. We see 
that for a 1.4M©, 1.4M© binary the minimum detectable 
values of A<& are ~ 8.1, 2.9 and 1.8 radians, for resonant 
frequencies of 16, 32 and 64 Hz. 

We conclude that it is unlikely that mode resonances 
will be detectable by LIGO, unless the spins of the neu- 
tron stars are anomalously large. Our results also show 
that resonances at higher frequencies are easier to detect. 



II. DETAILS OF ANALYSIS 

We start by reviewing the Fisher matrix formalism for 
computing parameter measurement accuracies; see for 
example Ref. [lj]. The inner product used on the vec- 
tor space of signals h(t) is given by: 



{hi I h 2 



k(f)h2(f) + hi(f)h* 2 (f) 

Sn(f) 



df, (2.1) 



where hi and h 2 are the Fourier transform of two signals 
hi and h 2 respectively, and S n (f) is the one-sided noise 
spectral density. We assume the model of the advanced 
LIGO noise curve used in [l[ 



Sn(f) = 



So [(fo/f) 4 



f < 10 Hz, 
2 (1 + (/V/o 2 ))] />10Hz 

(2.2) 

where Sq is a constant whose value is unimportant for 
this paper, and fo = 70 Hz. Given a signal model h = 
h(t,9' 1 ), where 9 l are the parameters of the signal, the 
Fisher information matrix is defined by 



dh 



dh 

del 



(2.3) 



and the rms measurement error in the parameter 9 l is 
V((A^) 2 ) = V^, (2.4) 

where X = T -1 . 

We next discuss our assumed form of the gravitational 
wave signal. In the absence of resonances, we use the 
model of Ref. 0]: 



h(f)=Af- r/6 e i *°W, 



(2.5) 



where 



$o(/3 = 2irft c -cb c -ir/4+-(8TrMf)- 5/3 
20 /743 11/A 



(2.6) 



Here M is the total mass of the binary, fi is the re- 
duced mass, M. is the chirp mass, and x = (nMf) 1 / 3 . 
We assume the signal shuts off at a frequency of / = 
(6 3 / 2 7tM) -1 . To include the effect of resonances, we re- 
place the phase $o(/) m Eq. (|2.5[) with the phase $(/) 
given by Eq. (|l.ip , which depends in addition on the 
two parameters fo and A<I>. [We approximate the band- 
width A/ rcs of the resonance to be zero.] Using this sig- 
nal model given by the seven parameters A, t c , (f> c , /i, 
M, fo and A$, we numerically compute the Fisher ma- 
trix (|2.3[) . and invert the matrix to obtain the parameter 
measurement accuracies shown in Table ILTl 

In these computations, we use the value A$ = 1. 
The phase perturbation will be detectable when A<I> > 
A(A$). From the form of the gravitational wave signal 



3 



it follows that the rms error A(A<f>) is independent of necessary for detection is given by the computed value of 
the value of A<£>. Therefore the minimum value of A<& A(A<£>), i.e. the last column in Table ITT1 
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